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SECTION  I 


INTRODUCTION 


In  many  applications,  thermal  energy  is  available  in  an  intermit- 
tant  form,  and  demand  for  this  energy  does  not  necessarily  coincide 
with  availability.  Solar  energy,  for  example,  is  not  available  at  night, 
and  similarly  wind  powered  electrical  systems  function  only  when  the 
wind  velocity  is  adequate.  Also,  the  greater  cost  of  commercial  electri¬ 
cal  power  during  periods  of  peak  demand  is  another  application  in  which 
energy  storage  could  be  used  to  good  advantage.  For  these  reasons, 
systems  which  can  economically  store  large  quantities  of  thermal  energy 
with  suitable  recovery  on  demand  are  of  increasing  interest. 

There  are  numerous  energy  storage  systems  available,  at  present, 
such  as  batteries,  hydrogen  storage/heat  engine,  flywheels,  etc.  Systems 
utilizing  the  phase  change  of  materials  having  a  high  latent  heat  of  fusion 
are  also  being  evaluated  at  present  and  are  of  particular  interest  because: 

1)  heat  transfer  occurs  at  nearly  constant  temperatures  and 

2)  large  quantities  of  thermal  energy  may  be  stored  in  a  relatively 
small  volume. 

The  basic  problem  with  systems  utilizing  the  phase  change  of 
materials  is  analytically  determining  the  melting  and  solidifying  char¬ 
acteristics  of  a  given  system.  The  problem  is  inherently  nonlinear  and 
other  than  a  few  exact  solutions  the  problem  must  be  solved  using 
numerical  or  analogue  techniques.  Carslaw  and  Jaeger^  discuss  the 

few  exact  solutions  available  as  well  as  some  approximate  solutions. 

(2) 

Kreith  and  Romie  present  solutions  for  the  cases  of  uni-dimensional 
cylinders,  spheres  and  semi-infinite  solids  for  the  condition  that  the 
temperature  gradient  in  the  solid  on  solidification  is  constant  or  that 
the  velocity  of  the  solid /liquid  interface  remains  constant.  Kreith  and 

l 


(2) 

Romie  also  discuss  analogue  computer  solutions  for  semi-infinite 
solids.  The  bulk  of  the  problems  involving  phase  change  have  been 
solved  by  numerical  techniques  and  include  a  variety  of  applications  and 

(3) 

problems  such  as  continuous  casting  of  cylindrical  ingots  ,  two- 

(4) 

dimensional  freezing  including  convection  effects  in  the  liquid  region  , 

and  axisymmetric  semi-infinite  hollow- cylinder  with  convective  heat 

(5) 

transfer  at  the  inner  radius 


SECTION  II 


OBJECTIVE  AND  SCOPE 

The  purpose  of  this  investigation  is  to  determine  the  thermal 
characteristics  of  energy- storage  systems  utilizing  a  cylindrical  annulus 
to  contain  a  phase  change  material.  The  outer  surface  is  insulated  and 
the  inner  surface  (a  heat  pipe)  is  subjected  to  a  constant  heat  flux. 
Additionally  the  influence  of  longitudinal  fins  within  the  annulus,  with 
regard  to  heat  transfer  characteristics,  is  to  be  assessed.  The  fins 
will  be  symmetrically  positioned  in  the  circumferential  direction. 

Since  a  solution  to  a  problem  in  solidification  is  also  a  solution  to 
the  corresponding  problem  in  melting  only  solidification  will  be  considered. 
The  system  is  of  sufficient  length,  in  conjunction  with  heat  pipe  charac¬ 
teristics,  that  axial  gradients  will  be  neglected.  The  phase  change 
material  will  be  assumed  to  have  a  distinct  fusion  temperature  and 
initially  the  liquid  is  uniformly  at  the  fusion  temperature.  Additionally 
the  thermal  properties  will  be  assumed  to  be  uniform  within  each  phase 
and  mean  values  will  be  used  for  each  phase.  Convection  within  the 
fluid  as  well  as  radiation  within  both  phases  will  be  neglected. 
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SECTION  III 


SYSTEM  DESCRIPTION 


The  thermal  storage  system  consists  of  a  cannister  surrounding 
a  heat  pipe.  The  annulus  formed  by  the  cannister  and  the  heat  pipe  con¬ 
tains  4.  1705  lb.  of  the  thermal  storage  material  which  is  LiF  -  MgF^  - 
KF.  The  cannister  and  heat  pipe  walls  are  made  of  .  065"  thick  304 
Stainless  Steel.  The  outside  diameters  are  1.  5"  for  the  heat  pipe  and 
3.  0"  for  the  cannister.  The  geometry  of  this  configuration  is  illustrated 
in  Figure  1.  Thermal  properties  assumed  for  the  materials  are  pre¬ 
sented  in  Table  I.  In  addition  there  are  two  finned  configurations 
containing  three  and  six  fins  of  1/16"  thick  304  Stainless  Steel.  The  fins 
are  parallel  to,  and  project  radially  outward  from,  the  longitudinal  axis 
of  the  heat  pipe.  The  configurations  with  three  and  six  fins  contain 
3.  9985  lb.  and  3.  8742  lb.  of  the  thermal  storage  material,  respectively. 
The  finned  configurations  are  illustrated  in  Figure  2. 
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TABLE  1 


HP/TES  MATERIAL  PROPERTIES 


Property 

Thermal 
Storage  Mat'l 
(LiF-  MgF^-KF) 

Fin,  Cannister  and 
Heat  Pipe 
(304  SS) 

3 

Density,  lbm/ft 

181.  0 

501.  12 

Conductivity,  BTU/ft  lb°  F 

4.  11 

13.  5 

2 

Diffusivity,  ft  /hr 

0.  035 

0.  224 

Latent  heat  of  fusion,  BTU/lbm 

350.  0 

N.  A. 

Fusion  temp,  F 

1310. 0 

N.  A. 

SECTION  IV 


ANALYTICAL  APPROACH 

In  order  to  evaluate  experimental  data  for  the  system  described 
previously  and  to  assist  in  the  design  of  similar  systems  it  is  desirable 
to  have  a  model  of  the  system.  The  information  sought  from  the  mathe¬ 
matical  model  is  the  position  of  the  liquid /solid  interface  and  the  tem¬ 
perature  distribution  within  the  solid  material.  Also,  the  time  of  final 
solidification  of  the  salt  is  sought.  For  the  axisymmetric  system,  the 
equations  describing  the  temperature  distribution  and  the  interface 
radius  are  solved  approximately  in  analytical  form  and  more  exactly  in 
numerical  form  by  solving  a  tridiagonal  matrix  eauation.  In  the  finned 
system  a  more  complex  formulation  involving  both  the  temperature  and 
enthalpy  is  used  because  of  the  complexity  of  finding  the  interface 
location. 
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SECTION  V 


ANALYSIS  OF  AXISYMMETRIC  SYSTEM 


Formulation  of  Equations 

In  the  axisymmetric  model  it  is  assumed  that  variations  along  the 
longitudinal  axis  of  the  system  are  negligible.  Heat  losses  from  the 
outer  shell  of  the  cannister  are  assumed  negligible  as  is  heat  transfer  by 
radiation  between  the  hot  outer  portion  of  the  system  and  the  cold  solidi¬ 
fied  salt.  The  entire  system  is  assumed  to  be  initially  at  the  fusion 
temperature  of  the  salt.  The  heat  flux  rate  at  the  inner  surface  is 
assumed  to  be  constant  with  respect  to  time  and  over  the  surface  of  the 
heat  pipe.  The  change  in  internal  energy  of  the  heat  pipe  wall  as  heat  is 
extracted  from  the  system  is  ignored. 


Under  these  assumptions  the  temperature  distribution  within  the 
salt  is  governed  by  the  unidimensional  transient  heat  conduction  equation: 


9T  _  5ZT  J_  ST 

at_Qf.  2'  RdR 
o  K 

where  T  =  temperature,  F 
t  =  time,  hr 
R  =  radial  distance,  ft 

2 

a  =  thermal  diffusivity  of  salt,  ft  /hr 


(1) 


Equating  the  heat  flux  rate  at  the  inner  surface  with  the  rate  at  which 
heat  is  conducted  into  a  unit  area  of  the  surface  yields: 


k 


ST 

SR 


for  R  =  R  ,  t  >  0 


(2) 


where  k  =  thermal  conductivity  of  the  salt,  BTU/ft  -  hr  -  F 
q"  =  heat  flux  rate,  BTU/ftZ  -  hr 
Rj  =  radius  of  outer  surface  of  heat  pipe,  ft 
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By  assumption  the  initial  temperature  distribution  is  described  by 


T(R,  t)  =  T  for  R  >  R.  t  =  0 

F  I 


where  T  =  fusion  temperature  of  salt,  F. 
r 

Also,  at  the  interface,  the  temperature  is  always  at  the  fusion  temperature 


T(R,t)  =  T  for  R  =  L,  t  >  0 

F 


(4) 


where  L  =  radius  of  liquid  solid  interface,  ft. 

Equations  (1)  through  (4)  constitute  a  partial  differential  equation  of 
second  order  in  space  and  first  order  in  time  along  with  the  necessary 
initial  condition  and  boundary  conditions.  However,  since  L  is  an  un¬ 
known  function  of  time,  the  problem  is  incomplete. 

To  determine  L,  two  approaches  are  considered.  In  the  first,  the 
heat  created  by  fusion  at  the  interface  is  equated  with  the  heat  conducted 
from  the  interface  to  yield: 

2irLphS£  ftT  =  2nL  k  for  R  =  L,  t  >  0  (5) 

3 

where  p  =  density  of  salt,  lbm/ft 

hg^  =  latent  heat  of  fusion  of  salt,  BTU/lbm 

Equation  (5)  along  with  the  initial  condition: 

L(t)  =  Rj  for  t  =  0  (6) 

is  sufficient  to  complete  the  formulation  of  the  problem  partially  posed 
by  Equations  (1)  through  (4). 

A  second  approach  to  determine  L  may  be  taken  by  equating  the 

change  in  internal  energy  of  the  salt  to  the  heat  conducted  out  of  the  salt 

at  the  inner  surface.  This  may  be  expressed  as 

R  =  L 

2  2  r 

PTT(L  -R  )hs  -  J  pc(T-T  \2TTRdR  =  q"t(2TTR  )  (7) 


where  c  is  the  specific  heat  of  the  salt. 


BTU/lbm°F. 


In  Equation  (7)  the  first  term  is  the  internal  energy  from  fusion  of  the 
salt,  the  second  term  is  the  change  in  internal  energy  from  the  solid 
salt  at  the  fusion  temperature  and  the  term  on  the  right  hand  side  is  the 
heat  extracted  from  the  inner  surface. 


There  will  be  a  temperature  drop  across  the  heat  pipe  wall  since 
heat  is  conducted  across  it.  Since  the  change  in  internal  energy  of  the 
wall  is  neglected  this  may  be  found  from  the  steady  state  heat  conduction 
equation  to  be: 


(8) 


where  AT=  temperature  drop  across  heat  pipe  wall,  F 

k^  =  thermal  conductivity  of  heat  pipe  wall,  BTU/ft  -  hr  -  F 
p  =  thickness  of  heat  pipe  wall,  ft. 

Systems  such  as  that  described  are  seldom  completely  filled  with 
solid  salt  in  order  to  avoid  excessive  hydrostatic  pressure  since  the 
salts  are  more  dense  in  solid  form  than  liquid.  The  ratio  of  the  volume 
of  solidified  salt  to  the  cannister  volume  may  be  written: 


(9) 


where  f  =  volume  ratio  of  solid  salt  to  cannister 

R  =  inner  radius  of  outer  cannister  wall,  ft 
o 

L,^  =  radius  of  liquid- solid  interface  when  solidification  is  complete. 
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To  generalize  the  solutions,  it  is  desirable  to  express  the  preceeding 
equations  in  terms  of  dimensionless  variables.  Substituting  the  variables 


0  =  k(T-TF)/q"Ri;  r  -  R/R^  i  =  L/R^  T  =  at/R2  (10) 

into  Equations  (1)  thru  (7)  yields,  respectively 


5  + 

J_  60 
r  6  r 

69 

6T 

(ID 

69 

6r 

1. 

for  r  =  1,  T  >  0 

(12) 

9(r,  T) 

=  0. 

for  r  >  1 ,  t  = 

0 

(13) 

9(r, t) 

=  0. 

for  r  =  l,  T  > 

0 

(14) 

M 

6t 

s  30 

STe  6r 

for  r  -  l,  t  > 

0 

(15) 

t(T)  = 

1. 

T  =  0 

(16) 

t2  - 
2  St 

1_ 

mt 

a 

j  9rdr  =  t 

(17) 

1 


Equation  (16)  may  be  solved  for  SL  to  yield 


1 


In  Equations  (15),  (17),  and  (18): 

sTe  =  q"RI/<*PhSje  (19) 

where  St  =  Stefan  Number 
Ae 

T  =  dimensionless  time 

=  dimensionless  frontal  radius 
r  =  dimensionless  radius 

9  =  dimensionless  temperature 


It  may  be  seen  from  Equations  (11)  through  (16)  that  the  solution 
for  6  and  &  in  terms  of  r  and  t  is  characterized  by  the  single  dimension¬ 
less  parameter  S^.  This  is  also  true  if  Equation  (18)  is  used  rather 
than  Equations  (15)  and  (16)  in  conjunction  with  (11)  through  (14)  to  for¬ 
mulate  the  problem. 


To  find  the  time  when  solidification  is  complete,  from  Equation  (9) 


and  from  Equation  (17) 


where  =  dimensionless  radius  when  solidification  is  complete 
T  =  dimensionless  time  when  solidification  is  complete. 


(20) 


(21) 


Solution 

Two  different  approaches  to  solving  the  system  of  equations 
representing  the  axisymmetric  system  will  be  considered.  In  the  first, 
the  change  in  internal  energy  of  the  solidified  salt  is  neglected  and  an 
approximate  analytical  solution  is  found  by  solving  Equations  (11)  through 
(15)  subject  to  this  assumption.  To  obtain  the  second  solution.  Equations 
(11)  through  (14)  are  expressed  in  finite  difference  form.  The  resulting 
tridiagonal  matrix  system  is  solved  in  conjunction  with  Equation  (18)  to 
find  the  temperature  distribution  and  interface  motion. 

Approximate  Analytical  Solution 

If  the  change  in  internal  energy  of  the  solid  salt  is  neglected 
in  comparison  with  the  heat  of  fusion  Equation  (11)  may  be  expressed: 
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I  ii 

r  Sr 


0 


(22) 


a2e 

Sr2 


+  ■ 


Then  the  expression  9  =  C  ^  In  r  +  is  seen  to  satisfy  Equation  (22). 
Applying  the  boundary  conditions  (12)  and  (14)  gives  the  solution: 

9  =  ln  (l)  (23) 

Substituting  (23)  into  (15)  yields 

s  a  _ 

St  "  ~  (24) 

In  Equation  (24)  we  may  separate  variables  and  integrate  to 
find  an  expression  for  JL; 


=  V 


1+2  St  T 


(25) 


Note  that  the  expression  (25)  satisfies  Equation  (18)  if: 


Aj 

/ 


9rdr  =  0 


(26) 


Equation  (26)  is  merely  a  restatement  of  the  assumption  that 
the  change  in  internal  energy  of  the  salt  is  negligible  which  will  later  be 
shown  to  be  quite  accurate  for  the  system  being  considered. 

Numerical  Solution 

The  analytical  solution  derived  above  for  the  axisymmetric 
system  is  an  accurate  solution  where  the  heat  of  fusion  is  numerically 
much  larger  than  the  specific  heat  of  the  storage  material.  If,  however, 
the  latent  heat  is  not  significantly  larger  numerically  than  the  specific 
heat,  then  the  assumption  made  to  obtain  that  solution  may  be  too  strong 
to  yield  realistic  results.  For  this  reason  a  numerical  solution  in  which 
the  internal  energy  of  the  salt  was  accounted  for  was  derived  and  is  pre¬ 
sented  here. 
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Let  the  temperature  distribution  at  time  t  be  known  in  the 

form  of 

k  k 

9i=0(r.,t);  l^i^n  (2 

kJ  1  k+ 1 

Then  if  ,  =  -£(t  )  is  known  we  may  find  9i  l^i^n+1  from  the 

k+ 1 

tridiagonal  matrix  equation  derived  in  Appendix  A  which  is  of  the  form 

'ai  bi  1  ‘9i  1  V 

C2  a2  b2  9 2  d2 

C3  *3  b3 


c  a 
n  n 


The  procedure  to  calculate  the  temperature  distribution  is 
illustrated  by  the  flow  chart  in  Figure  3.  First,  a  time  increment  At  is 
selected  and  r(2)  is  calculated  from  Equation  (25)  as: 


r ( 2 )  =  =  £(At)  =  t  1  +2St  At 

1  e 


6jl  =  S(r  =  l )  =  9(r(l))  =  2n(-|L) 


9  2  ■  ° 


where  the  9's  are  found  from  Equation  (23). 

k+ 1  k 

From  this  point  on,  the  procedure  to  find  9 i  =  9i  is  the  same  for 

k+ 1  k 

each  time  step.  Initially,  it  is  assumed  that  9i  =  0i  and  X,  is 

k+1 

calculated  from  Equation  (18)  where 


(32) 


/ 


9r<lr 


11 


i  =  2 


,  „k4  1  r 

U®:  'VVl 


-1 


>/2.  } ( r 


i-1 


)  +  ®_ 


kil 


n  !  1 


)/2. 


k4  1 

then  0.  ,  i=l,  n  are  calculated  from  Equations  (28).  This  temperature 

distribution  is  used  to  calculate  a  new  4  which  in  turn  is  used  to  calculate 
a  new  temperature  distribution.  This  iterative  procedure  is  continued 
until  the  change  in  4  meets  an  acceptable  error  criteria  at  which  point 
the  temperature  distribution  for  the  mew  time  has  been  calculated,  and 
the  procedure  is  initiated  for  the  next  time.  The  time  increment  is  held 
uniform  except  that  for  the  final  time  step,  is r  is  shortened  so  that  4(T^)  =  4^ 
from  Equation  (20).  The  computer  program  developed  to  calculate  the 
axisymmetric  temperature  distribution  is  documented  in  Appendix  B. 
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SECTION  VI 

ANALYSIS  OF  FINNED  SYSTEM 


Formulation  of  Equations 

An  approach  similar  to  that  used  for  the  axisymmetric  case  could 
be  used  to  model  the  finned  system.  However,  the  boundary  would  be 
much  more  difficult  to  locate  since  its  location  varies  with  the  angle  from 
the  fin  centerline.  Instead  an  approach  involving  both  the  enthalpy  and 
the  temperature  of  the  salt  employed  by  Shamsunder  and  Sparrow^  is 
used.  In  this  approach,  the  solid  liquid  boundary  is  found  as  one  of  the 
results  of  the  solution. 

Since  there  was  some  uncertainty  as  to  the  accuracy  of  this  method 
the  method  was  first  applied  to  the  axisymmetric  case  and  the  results 
were  used  to  examine  the  accuracy  of  the  method.  This  formulation  will 
be  presented  as  an  aid  to  understanding  the  approach.  Then  the  method 
will  be  used  to  formulate  the  finned  problem. 

In  this  approach  the  region  of  interest  is  again  divided  into  elements 
of  finite  size.  If  there  are  no  sources  of  energy  inside  the  control  volume, 
pressure  does  not  vary  with  time,  and  no  external  work  is  done  on  the 
control  volume  the  net  rate  of  increase  of  internal  energy  must  equal  the 
rate  at  which  heat  is  conducted  into  the  control  volume.  Then  enthalpy 
and  temperature  are  related  in  the  law  of  conservation  of  energy  for  the 
element. 


_d_ 

dt 


k  grad  T  •  ndA 


(33) 


where  h  =  specific  enthalpy  BTU/lbm 

n  =  unit  vector  normal  to  the  surface  and  directed  outward  from 
the  control  volume. 
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This  may  also  be  written  in  terms  of  dimensionless  variables 


71 


1 

p  Av 
s 


P(h-hx  ) 
hs4 


dv 


dimensionless  enthalpy 


(34) 


kf(T-T*) 

9  =  - : — t-  1  dimensionless  temperature 

Q  p  hs i  r 

f  s 

and  enthalpy  may  be  related  to  temperature  by 
fc(T-T-p)  ,  T  <  Tf 

*  I 

h-h 

[cCT-T^+hsX  ,  T  >  Tp 
Axisymmetric  Formulation  and  Solution 

Let  the  region  <  R  <  be  divided  into  n  elements,  each  a 

concentric  ring  of  width  AR.  Let  the  thickness  of  the  cannister  wall  be 
n^AR.  Then,  initially  T  =  Tp.  and  the  thermal  storage  material  is  all 
liquid  so  that  initially 

9=0.  i  =  1,  m 

Ti  =  0.  i  £  n  (37) 

71=1.  n^.<i£m 

where  m  =  n  +  n  total  number  of  nodes. 

Writing  (33)  in  terms  of  the  elements  gives 


(35) 


(36) 


St 


dt  p  2n|IAn 


/ 


«•? 


dn 

h*» 


i.  -  ttUwn-Vi 

I  Ar  2  2  ir  2 


(k  .  k  »(it.ir-TU/»H..<  .n|n 

I  AR  22  AR  2k, 


(38) 
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where  kj  =  thermal  conductivity  at  outer  edge  of  element 
k^  -  thermal  conductivity  at  inner  edge  of  element. 

Expressing  (38)  in  terms  of  the  dimensionless  variables  r,  t,  6,  and  11 
yields 

flr  Ml  =  (!+-*£_,  Ml  H-  — i  M  —  I 

dT  2ri  ks  I r  =  r.  ’  "  2r.  ’  ^ 

1  2  i  2 


which  can  be  expressed  in  terms  of  finite  differences 


-  <U2^'-2  <•*,> 

i  s 

If  the  element  is  considered  to  be  at  a  uniform  temperature,  substitute 
Equation  (36)  into  (34)  to  find  the  relation  between  temperature  and 
enthalpy: 


pc(T-T*) 


V  2 


.  ^  k  (T-T  )  a  k 

pc  s  s  .  s  „ 

Av  k  a  p  h«,  ark 

s  s  s  s 


0  =— 2-  T) 

or  k 
s 


T]  <  0 


=  0  T)  >  0 

Now,  rearranging  Equation  (39)  yields 


^  ^r-1  +  ,I+IT>  IT  9ui  ♦  T  #f.i  5  “2 - -  (41) 

■  ®  is 
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where  k,  =  <  2  k  k  /(k.4k  ) 
1  I  f  s  f  s 


k  =  <  2  k  k  /(k  +k  ) 
2  1  f  s  f  s 


h  =hf 
h  >h. 


h  sh 
h  =  hf  +  1 
h  >  h  +1 


Substitute  for  9  from  Equation  (40)  into  the  left  hand  side  of  Equation  (41) 
to  obtain 


[W  «  1!jl  {,14^  ^  *  d-|r>  r 3]  ^ 

L  at  “gk  2r.  ks  2r.  kg  J. 


Tlk  <  0 


Tlk  >  0 


At  the  inner  surface  of  the  innermost  node,  the  heat  flux  is 
known  so  that  Equation  (33)  may  be  written 


p  2nRARhs,a 


_ !LL*  m  =  _i  ii  ZTT(R+^)  «  P  hs£  -  2nq"R 

2  dT  k  SR  2  s  s  1  I 


where  q" 


(Ar)  /Tlk  _k-l  1  ,,  Arw.k  .  k  q"  ‘'I 

ljr(V\  1  =  -<1+2t:,<62-V  -  ^Thi:  4r  (43» 

c  1  co  “  y, 


s  s 


Notice  that  the  parameter  S^  has  again  appeared  on  the 
right  hand  side  of  Equation  (43).  Since,  in  the  previous  analysis  we 
concluded  that  this  was  the  driving  parameter  for  the  solution,  it  is 
logical  it  should  also  appear  in  this  formulation.  Equation  (4:>)  may  be 
rearranged  similar  to  Equation  (42)  to  yield 
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=  (4r£ 

At 


,  k-1 


k  .  ak 

IT 


kar  9, 


i+1 


-  St  Ar 


If  Tl,  <  0 


(44) 


Equations  (42)  and  (44)  are  solved  by  Gauss-Seidel  iteration 

with  over-relaxation.  The  procedure  is  to  first  calculate  the  right  hand 

k 

side.  Since  the  coefficient  of  T|.  is  always  positive  the  sign  of  the  right 
hand  side  is  the  same  as  that  of  Tj  so  that  the  appropriate  coefficient  may 
be  determined.  i  is  then  calculated  by  dividing  the  right  hand  side  by  the 
coefficient  of  "H .  When  the  temperature  distribution  has  been  found  with 
sufficient  accuracy  the  volume  solidified  and  the  change  in  internal 
energy  are  calculated  from: 


m 

vs  -  £  2rrr ^(T|°_max(Tl!<,  0.  ))  (45) 

i=  1 

m  u 

E,  =  E  2nr.(Tl°-D't)  (46) 

i=  1 

Finned  System  Formulation  and  Solution 

The  model  for  the  finned  system  is  set  up  very  similarly  to 
that  for  the  axisymmetric  system.  Since  the  temperature  may  vary  with 
angle  as  well  as  radius  two  spatial  dimensions  must  be  included  in  the 
model.  Because  of  the  symmetry  assumed,  only  the  shaded  region  in 
Figure  2  need  be  modeled  and  the  heat  flux  across  the  boundaries  =  0 

and  ♦  =  "/(no.  of  fins)  is  known  to  be  zero.  Also  because  heat  losses  are 

neglected  there  is  no  heat  flux  across  the  outer  surface.  The  heat  flux 
across  the  inner  boundary  is  specified  as  part  of  the  problem  so  that  the 
heat  flux  is  known  over  the  entire  boundary  of  the  region. 

We  write  the  enthalpy  equation  for  an  element  containing  a 

volume  fAv  of  fin  material  and  a  volume  (l-f)Av  of  salt  and  having  thermal 

conductivities  k  ,  k  ,  k,,  and  k,  on  the  surfaces  r  =  r  +  — .  *  =  *  + 

1  2  3  4  i  2  i  2  ’ 
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Ar  Aty 

r  =  r.  -  —  ,  and  <|i  =  t .  -  —  respectively.  Then  Equation  (33)  may  be 
written 

I*1-/  i,k  ,  III  ,  l2-  I  »JL  I*  ,  (i.— i  ^  .  -t  I  «.k 

At  i,j  I  ?.r  k  k  |rA*J  Ir  k  k  |  rA»  J  J  ifj 


lAlL  T)k-'  ,  {u*'-)  -±  ek  .  ♦  -i  -AL  0k 

AT  i,j  '  k  ill.)  k  rA"  I  •  J-*  I 


l-^r  e\ 


Now  we  use  Equation  (34)  to  establish  the  relation  between  T) 
and  9  assuming  the  element  at  uniform  temperature 


■fc  /// 


P(h-hs) 


pc  (T-T  )  p  c  T-T 

-f / r— 5  4  -L-5  -r-^-  (l-£)4v 

PgAvhg^  pg  Av  hs  ^ 


k  (T-T  '  )  a 

a  p  hSjL  ^  IT  PfCf  +  (1"f)^ 

s  s  s 


1  -  f(l--^-h) 
k  a  ' 
s  f 


T]  <  0 


0  <  T]  <  1. 


Equation  (47)  must  be  modified  for  boundary  elements  similar  to  the 
axisymmetric  case.  The  iterative  procedure  for  solution  is  illustrated 
by  the  flow  chart  in  Figure  4.  Here,  Equation  (47)  is  expressed: 

=  (DiT|i.y  4  c,,I’J)-T,Ul(j  ,c2,I'J>-7,u,i  +c3<I'J)-T,ik.I(j 
+  C4(I,  J)- Tlk  .  /C  5(  I.  J) 

',}  (4< 

C5(I,J)  =  D  MCI(I,  J)fC2(I,  JMC3(I,  J»+C4(I,  J))/t!-f  jl.j£pJ  ) 


when  T) . 


=  0. 
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INPUT 

INITIALIZE 


COMPUTE  M.N.D  IS  =  0 

COMPUTE  F,  1=  1 ,  M  J=  1 ,  N 
COMPUTE  K1 ,  K2,  K3 ,  K4, 1=  1 ,  M  J  =  1,N 
COMPUTE  C1,C2,C3,C4,C5,I=1,M  J=1,N 


1=1  J=1  ERR  =  0 


oh=i;1'  . 

>.  j 

COMPUTE  .  FROM  1491 


^ERR<EPS 


PRINT  IF  NEEDED 
k- 1  k 

1\  .  =  ~.  .  1=1.  M  j  =  l,N 

i.J  i.J 


Figure  4.  Flow  Chart  for  Enthalpy  Model  of  Finned  HP/TES  System 
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SECTION  VII 


RESULTS  AND  DISCUSSION 


Results  with  Axisymmetric  Case 

The  approximate  analytical  solution  as  previously  derived  in  dimen¬ 
sionless  form  is 


0  =  In  (r/X) 

i  -  Vi +  2  St  ~ 

e 

If  this  solution  is  correct,  the  change  in  energy  from  the  heat 
fusion  is 

2  2 

ttR  (i  -l)Ph8/  =  A  E. 

and  the  change  in  internal  energy  of  the  salt,  which  was  neglected  is 

Ro 


(23) 

(25) 


r  u  k 

AE2  =  /  ~ (T-Tp)  2nRdR 

R. 


-I 


1  2nq"Ri3  k(T-Tp) 


Of  q"R 


rdr 


I 


2nq"R 


3  l 


1  KI  r 

— -  /  0rdr 


(50 ; 


Substituting  Equations  (23)  into  (50)  and  integrating  yields 


2TTq"R3£2 


AE, 


'Rt£  r ,  ,2  T 


(51 
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The  ratio  of  the  neglected  energy  change  AE  to  the  change  from  the  heat 


fusion  AEj  is 


AE2  q"Ri  £2 


AE1  aphsA  A2  -  1 


Ln(A) 


=  ST 


e  A2  -  1 


Ln(  A) 


(52) 


In  the  worst  case  (heat  rate  =1.5  kw,  A  =  Af)  the  neglected  heat  is 
only  10.8%  of  the  heat  included  so  that  the  approximate  solution  is  fairly 
accurate  for  the  system.  Since  the  change  in  internal  energy  of  the  solid 
salt  is  such  a  small  part  of  the  heat  extracted  from  the  system,  the  soli¬ 
dification  times  predicted  by  this  simple  model  will  be  very  accurate, 
even  if  part  of  the  solidification  takes  place  on  the  outer  wall  due  to  heat 
losses.  In  this  case  however,  the  temperature  drop  would  be  less  than 
that  predicted  since  the  heat  of  fusion  would  not  have  to  be  conducted 
across  as  thick  a  layer  of  solidified  salt. 

Before  meaningful  results  could  be  obtained  for  the  numerical 
solutions  for  the  axisymmetric  case,  it  was  necessary  to  determine 
convergence.  In  this  program  the  increments  in  both  radius  and  time 
are  determined  by  the  parameter  At  and  was  for  each  case  half  that  of 
the  preceding  case.  Once  the  temperatures  agreed  for  consecutive  time 
increments  to  tnree  significant  digits,  convergence  was  assumed  to  be 
adequate.  For  the  runs  with  other  heat  rates,  At  was  adjusted  to  main¬ 
tain  the  same  increment  in  radius. 

In  Figure  5  temperature  is  plotted  as  a  function  of  radius  for 
several  different  times  for  the  axisymmetric  system  with  a  heat  extrac¬ 
tion  rate  of  1  kw.  Both  the  analytical  and  the  numerical  solutions  are 
shown.  It  may  be  seen  that  at  any  given  time  the  numerical  solution 
predicts  a  slightly  smaller  interface  radius  than  does  the  analytical 
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Figure  5.  Temperature  Distribution  As  A  Function  of  Radial 
Position  in  HP/TES  System  ( Axisymmetr ical  Case) 


solution.  This  is  because  a  small  amount  of  the  heat  extracted  comes 
from  the  change  in  internal  energy  of  the  solid  salt.  Therefore  the  salt 
solidifies  more  slowly  than  predicted  by  the  analytical  model  to  yield  the 
heat  which  is  extracted  from  the  system.  Also,  the  temperature  gradient 
at  the  inner  radius  is  less  steep  than  predicted  by  the  analytical  model 
since  less  heat  is  conducted  through  the  salt. 

Figure  6  shows  the  wall  temperature  plotted  vs.  time  for  the  axi- 
symmetric  system  with  several  different  heat  extraction  rates.  Also 
shown  in  this  figure  are  lines  of  constant  radius  of  the  liquid/solid 
interface. 

To  establish  convergence  for  the  enthalpy  model  of  the  axisymme- 
tric  case,  it  was  necessary  to  find  suitable  values  for  the  computer 
parameters  ERR,  NR,  and  NT.  See  Appendix  C.  It  was  found  that  when 
ERR  =  .  00001,  NR  =  63,  and  NT  =  80,  the  solution  agreed  within  .  002 
of  solutions  with  ERR  =  .  000001,  NR  =  126,  and  NT  =  160.  In  Figure  7 
the  approximate  analytical  solution,  and  the  numerical  solutions  with  the 
temperature  formulation  and  enthalpy  formulation  are  compared  for  a 
particular  heat  extraction  rate  and  time.  Agreement  is  seen  to  be  excel¬ 
lent  between  the  temperature  formulation  and  the  enthalpy  formulation. 
Even  the  slight  disagreement  shown  may  be  at  least  partially  explained  by 
the  fact  that  the  cannister  wall  is  included  in  the  enthalpy  model  but  not 
in  the  temperature  model. 


Results  with  Fins 

Because  of  the  run  time  and  memory  requirements,  the  conver¬ 
gence  criteria  for  the  finned  model  was  based  on  results  with  the 
axisymmetric  model.  The  convergence  study  on  the  axisymmetric  model 
showed  that  NT  and  ERR  were  most  detrimental  to  accuracy.  See 
Appendix  D.  Therefore  the  finned  model  cases  were  run  with  computer 
parameters  NR  =  21,  NS  =  30,  NT  =  80,  and  ERR  =  .  00001. 
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Figure  7.  Temperature  Distribution  in  Axisymmetric  HP/TES  System 
Determined  Analytically  and  by  Two  Numerical  Techniques 


Figures  8,  10,  12,  14,  and  16  show  the  interface  radius  plotted 
versus  angle  from  the  fin  centerline  for  heat  extraction  rates  of  .  25,  .  50, 

.  75,  1.  0,  and  1.  5  kw  at  several  times.  They  show  that  a  significant 
amount  of  freezing  takes  place  along  the  fin  but  the  close  agreement  of 
the  interface  shape  for  the  two  configurations  indicates  a  similar  amount 
of  heat  is  extracted  by  the  fins  in  each  case.  Figures  9,  11,  13,  15,  and 
17  show  the  wall  temperature  profiles  for  the  same  cases.  Note  that  a 
significant  reduction  in  the  temperature  drop  is  obtained  with  fins.  This 
reduction  is  seen  to  increase  with  time  and  with  heat  extraction  rate.  It 
appears  from  the  figures  that  most  of  the  temperature  drop  obtained  with 
6  fins  is  also  obtained  with  three  fins. 

Figure  18  shows  the  wall  temperature  as  a  function  of  the  heat 

extracted  for  the  zero,  three  and  six  fin  configurations  and  a  1  kw  heat 

rate.  The  effectiveness  of  fins  as  a  device  for  reduction  of  total  system 

weight  was  assessed  in  the  following  manner.  The  maximum  allowable 

dimensionless  temperature  drop  was  arbritrarily  selected  to  be  0.4. 

From  Figure  18  the  amount  of  heat  extracted  was  as  shown  in  column  2 

of  Table  2.  The  interface  radius  for  the  axisymmetric  configuration  was 

^  =  1.513  and  with  a  30%  void  volume  allowance,  the  outer  radius  is 

r  =  1.6857  and  the  dimensionless  volume  is  v  =  1.8416.  The  other 
o 

cases  were  sized  to  provide  for  the  same  volume  of  salt  resulting  in  the 
outer  radii  shown  in  column  4  of  the  table.  The  lengths  required  to  pro¬ 
vide  an  equivalent  amount  of  heat  are  shown  in  column  3  and  the  dimen¬ 
sionless  system  weights  are  shown  in  column  5.  It  may  be  seen  that 
despite  the  shorter  lengths  for  the  finned  configurations,  the  systems  with 
fins  are  much  heavier  than  the  conventional  configuration  without  fins. 
Thus  the  given  fins  are  probably  not  an  effective  way  to  reduce  the  tem¬ 
perature  drop  when  system  weight  is  a  prime  consideration. 


HP/TES  Systems 


DIMENSION 


DIMENSIONLESS  WALL  TEMPERATURE,  «  =  K(TM-TF)/q 


ANGLE  FROM  FIN  %  ,  DEGREES 


Figure  13.  Wall  Temperature  vs.  Angle  for  0  and  6  ' 
HP/TES  Units,  Heat  Extraction  Rate  =  0. 


37 


DIMENSIONLESS  WALL  TEMPERAT 


Figure  15.  Wall  Temperature  vs.  Angle  for  0,  3,  &  6  Fin 
HP/TES  Units,  Heat  Extraction  Rate  =  1.00  kw 


DIMENSIONLESS  HEAT  EXTRACTED/UNIT  LENGTH,  Q/tt  R*  /> h 
0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 
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TABLE  2 


FIN  EFFECTIVENESS 


Heat  Extraction  Rate  -  1  k\v 
Dimensionless  Wall  Temperature  -  -0.4 


No.  of 
Fins 

Heat  Extracted 

Unit  Length 

Relative 

Length 

VRi 

Relative 

Weight 

0 

1.  36 

1.0 

1.686 

1 . 

3 

1.52 

.  8947 

1.766 

1. 0387 

6 

1. 62 

.  8395 

1. 823 

1. 1 189 
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SECTION  VIII 


SUMMARY  AND  CONCLUSIONS 

The  heat  transfer  process  has  been  successfully  modeled  for  both 
axisymmetric  and  finned  systems.  For  axisymmetric  systems  the  tem¬ 
perature  formulation  may  be  solved  very  quickly,  providing  a  solution 
which  is  extremely  accurate.  The  difficulty  in  locating  the  interface 
boundary  makes  it  difficult  to  adapt  this  formulation  to  two  dimensional 
systems. 

For  those  systems,  the  enthalpy  model  provides  a  system  in  which 
the  temperature  distribution  and  interface  boundary  are  found  simultan¬ 
eously.  When  this  formulation  was  applied  to  the  axisymmetric  system 
it  was  found  that  accuracies  were  as  good  as  those  attained  with  the 
temperature  model  but  that  an  order  of  magnitude  greater  computer 
time  was  required  to  obtain  the  solution. 

For  the  finned  system,  an  additional  order  of  magnitude  of  computer 
time  was  required  to  obtain  the  solution.  The  solutions  with  fins  are 
slightly  less  accurate  than  the  axisymmetric  since  it  wasn't  practical  to 
study  convergence  fully.  Solutions  were  attained  with  sufficient  accuracy, 
however,  to  give  a  meaningful  comparison  of  the  finned  systems  to  the 
axisymmetric  system.  It  was  found  that  some  reduction  of  the  tempera¬ 
ture  drop  at  the  inner  wall  was  attained  with  fins  but  that  systems 
employing  fins  would  in  turn  be  heavier.  Thus  fins  may  not  be  practical 
for  this  particular  fin/system  design. 
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SECTION  IX 
RECOMMENDATIONS 

Although  the  analytical  and  numerical  techniques  developed  were 
very  accurate,  the  final  solution  may  be  adversely  affected  by  the 
accuracy  of  the  measured  properties  of  the  salts  as  well  as  the  limiting 
assumptions  used  in  deriving  the  governing  equations.  Based  on  the 
analysis  completed  it  is  recommended  that 

1)  The  salt  physical  properties,  liquid  and  solid,  should  be 
accurately  determined. 

2)  The  effect  of  radiation  to  and  from  the  salt,  liquid  and  solid, 
should  be  assessed.  Radiation  properties  of  the  salt  are 
needed. 

3)  The  experimental  apparatus  should  be  redesigned  to  (a)  relieve 
hydrostatic  stresses  on  melting,  (b)  improve  calorimetry  and 
(c)  determine  a  fin  configuration  to  optimize  heat  transfer. 

4)  The  temperature  model  for  the  axisymmetric  case  may  be  use¬ 
ful  to  study  boundary  conditions  which  vary  with  time  since  it  is 
solved  so  quickly.  It  is  recommended  that  this  program  be 
suitably  modified  and  studies  of  the  effects  of  time  varying  heat 
extraction  rate  be  conducted. 

5)  The  enthalpy  model  has  the  versatility  to  study  solidification 
fronts  in  unusual  geometries.  This  model  should  be  refined  to 
shorten  the  solution  time  and  better  define  the  accuracy.  It 
could  then  be  used  to  examine  other  two  dimensional  configura¬ 
tions,  including  the  axial  variation  of  the  nonfinned  case  and  the 
assymmetry  due  to  the  volume  void  on  solidification. 

6)  The  model  of  the  axisymmetric  system  should  be  expanded  to 
include  the  heat  transfer  within  the  heat  pipe  and  the  calorimeter. 
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Such  an  expanded  model  would  be  useful  in  the  design  of  further 
experiments  with  such  systems. 
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APPENDIX  A 


AXISYMMETRIC  HEAT  CONDUCTION  EQUATION 

In  this  appendix  the  axisymmetric  1 -Dimensional  heat  conduction 
equation  in  cylindrical  coordinates  is  written  in  finite  difference  form 
for  a  variable  mesh  grid.  This  equation  is: 


s 9  j_  _s©  _  _s© 

.2  r  Sr  St 

Sr 


(A.  1) 


Where  0  is  the  dimensionless  temperature  at  time  t  and  radius  r.  The 
object  is  to  express  Equation  (A.  1)  at  r  =  r  in  terms  of  the  temperatures 

K 

9k_1'  9k  and  9k+1  at  radii  rk_  j  *  r  .  and  r  respectively.  To  do  this 
write  0^+1  and  by  Taylor  Series: 


^  9 

9k+i  =  +  (rk+rrk'57  + 


(rk+l"rk*  S28 


2  ! 


Sr 


e  =0 

k-1  k 


,  ,  se  ,  <rk'rk-i’  a2e 

(Vrk-l>S7  +  - 2! -  —  + 


dr 


(A. 2) 


2 

d  0  u  q 

Neglect  higher  order  terms  and  solve  Equations  (A.  2)  for  —  and - 

or  2 

to  obtain  "r 


ae  w  (Vrk-i>Z  9kki  4  (rkH  I'Vi1  9k  -  (rk,i  -k 

dr 


-r  J2  V 


,rk<rrk"Vrk-i,<rk+rrk-i» 


1^.  2  'VVi'  Vi  -  (rkirrk-il  ek 4  (rkki-y  9k.i 
ar2"  (rk  +  l-rk),Vrk-l>,rk+l-rk-l) 


(A.  3) 
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Substitute  Equation  (A.  3)  into  the  right  hand  side  of  (A.  1)  and 
group  terms  to  obtain: 


a  e  j_a_e  _  _ 1 _ 

a,2  *  r3r  "  rk|rk*rrl.i1('nrrkllV'k.il 
t<3Vrk+i>  ''k+r'k1  \-i +  'Wr'k.i1  (rk+rIk-i)  ei 


+  (3Vrk-I,,Vrk-l,ek+l) 


■  2(  c  0  -fa  0  +b  9 
'  k  k-1  k  k  k  k+lJ 


(A. 4) 


Let  9  represent  the  temperature  at  time  T  and  ♦  represent  the 
temperature  at  an  earlier  time  T-At.  The  right  hand  side  of  Equation 
I  A.  1)  is  written  by  finite  differences  and  the  Crank-Nicholson  method 
is  employed.  Then: 


1  fa20  +  i  s9_  ^  a 

2  S  2  r  dr  2  rdr 

Idr  Sr 


=  c  (9  +♦  ) 

ky  k-i  k-r 


9  -<t> 

+  a'  (9,  +♦,  )  +  b,  ( 9  +♦  )  =  =  -7  k 

k'  k  k  k  k-f  1  k+1  St  At 


(A. 5) 


Thus,  grouping  the  unknown  terms  on  the  left  hand  side  the  finite  differ¬ 
ence  form  of  Equation  (A.  1)  is 


c  9  ,  -f  a  0  +  b  9  =  d. 

kk-1  kk  kk+1  k 


(A.  6) 


where 


(3rk"W  (rk,rrk> 

2rk(rk+rrk-l)  <rk(  l"rk*  |rk'rk-l! 


(3Vrk-1)(rk-rk-l) 

2rk(rk+i'rk-i’  (rk+rrk)  (Vk.i' 


I 


'“Vk-i-W  <rk+rrk-i' 


^k^k+r'k-i’  (rk+rrkl  !rk'rk-i)  lT 


'k  k-1 


(-a,  4— ^ -)♦,  -  b  * 

k  At  k  k  k+1 


If  is  known  for  k  -  1 , 2,  .  .  .  , n  where  r  is  the  radius  of  the 
K  n 

interfacial  surface  at  time  t-At,  Equation  (A.  6)  is  applicable  for  k  = 

2,  3,  ....  n-  1 .  At  the  inne r  boundary  an  imaginary  node  is  ass umed  s uch 


r2  rl  rl  r0 


(A.  7) 


and  the  inner  boundary  condition  is  expressed  in  finite  difference  form: 


v*o  .  V!o 
2(Vl>  ‘  2(r2-rl> 


(A.  81 


substitution  of  Equations  (A.  7)  and  (A.  8)  into  Equation  (A.  6)  for  k  =  1 
yields  the  equation  for  the  inner  boundary: 


{-S  --i- 

I'W  J 


1  e2  =  — * — ?  +  ••  -  — —  «  — —  .A. 9) 

)  <r,-r,)Z  2  1  V.  Vrl 


For  k  -  n  Equation  (A.  4)  is  not  valid  for  $  since  - —  is  not  defined  at  r  . 

or  n 

Instead,  from  Equation  (A.  1) 

2 

5  ♦  ,  d$  d$  Sr  d$  _  , 

_ n_  +  _1  _ n  _  _ n  _  _ n  _ n  _  n  d^ 

-.2  r  dr  dT  dr  dT  dr  dT 
o  r 

and  from  Equation  (15) 


n  _  _ 1_  d / 

dr  2ST  dT 


so  that,  in  finite  difference  form 


49 


2 

a  $ 

n 

,  d* 

1  n  - 

1  i-r 

ar2 

1 

S' 

1 

iTn-fl 

(A. 7) 


and  substitution  into  Equation  (A.  5)  yields 


(2-r  )  * 

c  9  J-  a  9  =  - - - —  /2  S„  + 

n  n-i  n  n  _t  2  Te  At 

n+1  n 


(A.  8) 


'--ere  c  and  a  are  given  by  Equation  (A.  6)  with  l  used  for  r 

n  n  k+1 


50 


APPENDIX  B 


PROGRAM  ICE2 


This  appendix  documents  the  program  ICE2  which  finds  the  tem¬ 
perature  distribution  and  interface  radius  for  the  axisymmetric  system 
using  the  temperature  formulation  developed  in  Section  V.  2.  b.  Figure 
B-l  is  a  program  listing.  Figure  B-2  shows  the  necessary  control 
cards  and  sample  input  data  to  run  the  program.  The  input  is  organized 
as  follows: 


Record  1  has  a  format  of  5E12.0  and  contains  the  following  variables 

2 

DIF  -  diffusivity  of  thermal  storage  material,  ft  /hr 

CON  -  thermal  conductivity  of  thermal  storage  material,  BTU/ft-hr-°F 

3 

DEN  -  density  of  thermal  storage  material,  lbm/ft 

TF  -  fusion  temperature  of  thermal  storage  material,  °F 

HSL  -  latent  heat  of  fusion  of  thermal  storage  material,  BTU/lbm 

Record  2  has  a  format  of  5E12.0  and  contains  the  following  variables 


RI 

RO 

QR 

DNR 

DT 


-  outer  radius  of  heat  pipe  wall,  ft 

-  inner  radius  of  outer  canister  wall,  ft 


-  heat  flux  rate  at  outer  heat  pipe  wall,  BTU/ft  hr 

-  ratio  of  solid  salt  volume  to  canister  volume 

-  time  between  successive  calculations, 


'I 


Record  3  has  a  format  of  4E12.0  and  contains  the  following  variables 


2 

DFF  -  canister  wall  diffusivity,  ft  /hr 

CFN  -  canister  wall  conductivity,  BTU/ft-hr-°F 

3 

DFN  -  canister  wall  density,  lbm/ft 
P  -  canister  wall  thickness,  ft 
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Record  4  has  a  format  of  El 2.  0,  315  and  contains  the  following  variables 

EPS  -  tolerance  on  calculation  of  energy  integral 
IP  -  print  every  IPth  time  calculations  are  made 
L  -  if  L  =  1  output  is  dimensionless  otherwise  it  is  dimensional 
ICN  -  case  number  for  identification 
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APPENDIX  C 


PROGRAM  ENTH 


This  appendix  documents  the  program  ENTH  which  calculates  the 
axisymmetric  temperature  distribution  and  interface  radius  from  the 
enthalpy  formulation  developed  in  Section  VI.  l.a.  Figure  C-l  is  a  pro¬ 
gram  listing.  Figure  C-2  shows  the  control  cards  and  sample  input  data. 
The  input  is  organized  as  follows: 

Record  1,  format  5E12.0 

2 

DIF  -  thermal  diffusivity  of  thermal  storage  material,  ft  /hr 

CON  -  thermal  conductivity  of  thermal  storage  material,  BTU/ft-hr-°F 

3 

DEN  -  density  of  thermal  storage  material,  lbm/ft 

TF  -  fusion  temperature  of  thermal  storage  material,  °F 

HSL  -  latent  heat  of  fusion  of  thermal  storage  material,  BTU/lbm 

Record  2,  format  5E12.0 

RI  -  outer  radius  of  heat  pipe  wall,  ft 

RO  -  inner  radius  of  outer  canister  wall,  ft 

2 

QR  -  heat  flux  rate  at  outer  edge  of  heat  pipe,  BTU/ft  hr 

DNR  -  ratio  of  solid  thermal  storage  material  volume  to  canister  volume 

ERR  -  allowable  error  in  enthalpy  change 

Record  3,  format  4E12.0 

ORF  -  relaxation  factor  effects  speed  of  convergence 

2 

DFF  -  thermal  diffusivity  of  canister  wall,  ft  /hr 

CFN  -  conductivity  of  canister  wall,  BTU/ft-hr-°F 

3 

DFN  -  density  of  canister  wall,  lbm/ft 
Record  4,  format  715 

NR  -  number  of  nodes  in  radial  direction  in  salt 

NT  -  number  of  points  in  time 


IM  -  maximum  number  of  iterations  allowed 
IP  -  print  every  IPth  time 
IR  -  dummy  variable 

IF  -  number  of  nodes  in  radial  direction  in  wall 
ICN  -  case  number  for  identification 
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.  S ?- 1 ,  Figure  C-l.  Program  ENTH 
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IF  U.CE.IF)  af=uff 
if.  (I.LT.IFI  Dl.l  =Cf N/CCN 
If— 43  .CE.IF4-  6142-Cr^l/CON 

If  II.EO.IF)  CK1  =  ?.*CCN«  CFf,/  UCN+CfN) /CC'. 

-  -if  -<  t  »  f  Q-»  I F  ♦  1 )  - - 

H(;  )  =  CD*FH  <!)♦(!  .♦CII)  >H  M«T  <  1*1)  ♦  <  l.-c  (1  )  l-CN2*T  II-ll|/n*ICM 

-  - - - - 

IF  <H<I)  ,LT •  0. )  GO  TO  30 

hc:»=h«i  imdmcni*  <i.*c(i  n*ci.?*  n.-.iini/on/D 

- 7  <14=0. -  —  - 

GD  TO  40 

—  70  -  -7-3r-?26yE-  - -  -  - -  -  . 
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sW=ff.+  :  C?  (H(I  J-CH) 
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20  00.  r  IKl'E 
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APPENDIX  D 


PROGRAM  FIN 

This  appendix  documents  the  program  FIN  used  to  calculate  the 
temperature  distribution  and  interface  radius  for  the  finned  HP/TES 
systems.  Figure  D-l  is  a  program  listing  while  Figure  D-2  shows  the 
control  cards  and  sample  input  data.  The  input  data  is  organized  as 
follows : 

Record  1,  format  5E12.0 

2 

DIF  -  thermal  diffusivity  of  salt,  ft  /hr 

CON  -  conductivity  of  salt,  BTU/ft-hr-°F 

3 

DEN  -  density  of  salt,  lbm/ft 

TF  -  fusion  temperature  of  salt,  °F 

HSL  -  latent  heat  of  fusion  of  salt,  BTU/lbm 

Record  2,  format  5E12.  0 

2 

DFF  -  diffusivity  of  canister  wall,  ft  /hr 

CFN  -  conductivity  of  canister  wall,  BTU/ft-hr-°F 

3 

DFN  -  density  of  canister  wall,  lbm/ft 

2 

QR  -  heat  flux  rate  at  outer  edge  of  heat  pipe,  BTU/ft  hr 
DNR  -  ratio  of  solid  salt  volume  to  canister  volume 

Record  3,  format  5E12.0 

RI  -  outer  radius  of  heat  pipe  wall,  ft 

RO  -  inner  radius  of  outer  canister  wall,  ft 

B  -  half  thickness  of  fin,  ft 
P  -  thickness  of  heat  pipe  wall,  ft 
G  -  length  of  fins,  ft 

Record  4,  format  2E12.0 

ORF  -  relaxation  factor  effecting  speed  of  convergence 
ERR  -  error  tolerance  on  enthalpy  change 
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Record  5,  format  415 


NR  -  number  of  nodes  in  radial  direction 

NT  -  number  of  points  in  time  ;  i 

NS  -  number  of  nodes  in  angular  direction  ( 

NFN  -  number  of  fins 

Record  6,  format  415 

IM  -  maximum  number  of  iterations 

/ 

IP  -  print  every  IPth  time  calculated 

ID  -  print  diagnostic  if  ID  =  1 

ICN  -  case  number  for  identification 


*  I*S»  tSS*SSt*3?S3m?t*S?tiSSSt?SS?SS*?StSSStSSS?S?SSSS?*:tT.?.*STT?SSST.S 

'  1 1 :  S  J 1  i  ?  C  -  5  ^  J  i  »  1  3  f  ?  ;  i  ST  -  i  J  l  Z  Z  :  i  *  i  2  i  £  I S  j  ?  ;  _  i  *  :  Z  ;  S 1  ?  ?  i  i  ;  3  $  J  ?  a  ?  t  f I  l 


»M  fp. 


7*4/71, 


CP,=1 


F  T  U  4.6+  446 


CnXG*.*, 
ul  *£  KS 
1»H (25, 
-  £  t  r  c 
>E-r  ( 


r  rlN(I  »iPU  1  ,  T  APE  1  =  J  N  F  JT  ,  OU"  P  UT  ,  T  APE  2=OUTPU7  ) _ 

I  OV  C 1 <  25 , 30 ) , C  2  < 2  5 ,3 0  I  ,  C3  <  25,  30"  1  ,  C4  ( 25 , 3  0  )  ,  C  ~5  (  25,  30  )  , 
30) , T (25, 30) ,  H<25, 30) ,NF (25) ,F (26) ,RF (30) 

1,10  0)  OIF  ,C-N,C":S.  ,1f  ,PiL 

i.ioo)  c,rF  ,cr.N,cPt.,of  .pi.t 


<1,100) 

(1,100) 

(1,100) 

(1,100) 
(  1 , 20  0) 
( 1 ,200) 


C.rF  ,CF.-(,Lrl.,Of 
Ri  ,  f\0  ,  3  ,  P 
0-F 

■IP  , T  ,  *, S  , 
iK,4.p,ir , 


J?=l 
=0  =  5  7 

O'  =  (  . 

0*  -3.1  4156  27/^FN/N'3 
37  -Or.*  I  :/C1F/D£N/H3l 
Z  r  =  (  X 


I  P  <T  .LE.1F  )  GO  70  1 

;fi;.g*,k-if)  c-c  to  10 

u=  (2.*  (I-iF)-l .  )  / 2.-»0P>7fI 
.1=5/  PI  /4/or 


J.F  (iC  .Li  .  1)  GO  10  4 

W-  PL  (2,5  0  0^  f  ^-7  ,  ,  ET  ,  C  FP 

win  e  c  2 ,  ci  o )  _rr,  rs ,  1  ftm  ,  p  1 ,  *• ,  m 

Wfl’i  (2,700)  (F(1),I=1,K) 


WPPE  <2, £00)  ( N  F  <  1 )  »  I  =  1  »  K ) 


W°.1T  E  (2, £00) 
CONTINUE 


00  20  1=1, r 
d:  20  j=i,* 


Figure  D-l.  Program  FIN 


MIS  7AM  I»  B»T  WAA.ITI  ******** 
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Ut  tf* 


:ic .  j  >  =  c?  <: ,  j>=o3  <i ,  j  >  =  c«*  i:.  j  >  =  coh 

IF  (  J.LE.NF  (I)  )  Cl  (I ,  J)  =  G2  (I.  J)  =  C3(I  ,  J)  =  C4  (I ,  J)  =CFK 

IF  (  j.Fn.tFIDHI  l  1  ( I  *  J )  =TTf  ]  ,  J  r=  ( 1-T~  fTTT  ~*T.  JF+TTT  rrC  FT! - 

_ TF  <  J.  FO.t.F  (j  )  .  :;g.r  II  )  .IT  .  .5)  C2  <1,  J)  =  CC,'-'»CFN/  (  .5*  tC  C  M+CF  N)  +  F  (  I  ) 

1*  (  -.OIJ-CFM  )  ~  “ 

_ IF  (  J.EO.NFtH*l.ANO.F(I)  I  .  .5)  C4(I , J)  =C  ON*  C  F  H/ (  .  5*  (  COM  CFN )  ♦  F  (! 

1)* CCOh-CFK) )  "  “ 

_ IF  tJ.F0.NF(I)*l.&N0.F(I).GE..5)  G2(l,J)=t0N*CFN/(.5M3.*CFH-C0N)  * 

_ iF(n*c;ct.-cFN)) 

IF  (  J.FC.Nc  (I)  +1  .AN3.F  (I  >.GS.  .  5)  CMI,  J)=CFF> 

IF  (  J.EO.KF(I)*?.AND.F  <1  ).GF..5)  C  4(  I ,  J)  =CON*  CFN/ l  .  5*  (  3.  *C  FN-CON  )  ♦ 
IF  (I  )  *  (uC.-CF  J)  ) 

20  CCMIMJE 

IF  (ID. LT . 1)  GO  i C  6 

W^ITE  (2,700)  < (C1(I,J), J=1,N)  ,I=1,M) 

W~.  I"  E  (2.700)  (  (C2C,  J),  J=  1 » »l »  ,I=1,M> 

_ M~;7t  <  2.700)  (  (i.3(I,  J).  J=1,N)  ,1  =  1, H) _ 

W-.ITE  (2,700)  (  (C4(I,J),J=1,N>  *1  =  1, M) 

6  COM  IKUE 

DC'  3  0  I  =  1 ,  F. 

1  =  1.  M2.*  C-IF)-1.)/2.*DR 
DO  30  J=1 , K 
IF  (*»EQ«)*)  GO  TJ  40 

Cl U , J  )  =  2.* Cl (I , J)*C3 (1*1, J)/ (Cl <1,  J)*C3 (I ♦! , J) ) 

_ :3ii»i,  j)=cki,j) _ _ 

40  :i(I,J)=J  ( I  *  J ) /O  ON*  (1  .♦CF/2./F) 

0 2  ( I  ,  J  )  =C  2  ( I ,  J )  *  (  Dr,/(./DS  )  * '2  ./CON 
03  C  ,  J)=^3  (I,  J)/uOi<*  (1.-CF/2./K) 

:4C,J)=G4(I,J)*<CK/f /DS)**2./CCN 
IF  C.FO.l)  C  3  ( I ,  J)  =  0. 

IF  (J.EO.l)  C4(I , J)=0. 

IF  C.EO.K)  C1(I,J)=0. 

IF  (J.EQ.N)  C2C,J>=0. 

Fr  =  0. 

_ IF  (J.i.E.NF(I))  FF=1. _ 

IF  < J.EQ.KF(I>*1)  FF=F (I ) 

05(1,  J)=C+  (Cl  (I,  J)+02(I,  J)+C3C,  J)+C4(I,  J)  )  /  ( 1. -FF*  C 1 .  -DFP  ) ) 

~C,  J)  =  0. 

_ 3H(T  , J)  =  l. _ 

IF  (J.IE.NMIII  PH(I,J)=0. 

IF  (  J.'Q.nF«1)«’1)  PH(I  ,J)=1.-F  (I) 

MCI, J)=PH(2,J) 

30  OC-::  IMIE 

TF  ( I D  •  l .  1 )  Gu  *.  0  7 

i  (2,700)  (  LKIi  J)  ,C2  (I,  J),C3(1 ,  J)  ,C4(I  ,  J)  ,F5(I,  J)  , 

1J=1,N> ,1=1 .*) 

w-.i-E  (2,700)  iihc.ji  ,j=i,n) ,:=i,m> 

W  -  I T  E  (2,700)  ((FH(I,J), J=1,N) ,1  =  1, M) 

WrM7£  (2,700)  (  (TII,  J)  ,J=1,I.)  ,1=1, Ml 

7  COMISUE 
,N.-.=  0N«*100 

M--  I7E  <  2 , 3 00)1  Clt , KFN,  ,"Ur dT^TTB , p- 0 ,  G 

bi  -ONF/100 

WPITE  (  2,350)  DEf; ',OFN,C  ON, OF^i, OIF, OFF, HSL,H,K,fJT  ,OFF,rRF,If 
_ W-.ITE  (2,.3fe0)Tr.EF,*’F,ST  _ 
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10  6  J=l,6 
RF  <  J)  =  0  • 

0H=H(1,1) 

H  <  1 ,  1 )  =  <  3*  P  H  <  1 , 1 )  ♦  SI  ( 1  , 1  )  *  T  ( ?  ,  1  )  ♦ 
Ml 1,11  -0SF* ( H ( 1 « 1 ) -OH ) +0H 


T (1,1) =H(l,l)/CFn 
EN  =  /.6S(H<1,1)-CH) 


ED=m9S (Mil, II) 

DO  5  0  1=?,M1 _ _ 

OH=H(T,i) 

M ( I ,  l )  =  <D*  PH  <  1 , 1 )  *01  ( 1 , 1 )  *^TTi  1  ,TT*C 2  C ! ,  1)  *T  ( I ,  2T*T3n7l7* 

1) ) /C5( 1 , 1) 

IF  (H(I,1) .Lt.O. )  GO  TO  60  “ 

H(T  ,1)  =H  (T  ,  1)»C5  (1  ,1 )  /b 


IF  <r.F  ( 1 )  .LT.l.  )  RF(1)  =  (I-IF-H(I  »  J ) )  *  D  R*  1 
GO  TO  FO  ~ 

IF  (OH.LT.C.)  h(A,l)=0‘F*(H(I,l)-0H)*0H 
FF  =  F  ( I  ) 

IF  (NF(I).GT.O)  FF=1. 


» 1 ) =  H ( I , 1 ) / < 1 •  — FF*C 1. -DFR) ) 
■£')♦  '.FS  (HCI  ,1 )  -Om> 


l»=Eu*ASS(MI,1)  ) 

^  1  ItiUE 
'"  H=H  ( f'  ,  1 ) 

1  )  =  <0*F  H(M,  1)  +02  (K,l  )  *T  (“  ,2)  +C3<  K,  1)  *T  (M-l,  1  )  )  /C5  <M,  1  ) 


-S  *  SI )/C5(l,J> 

Mil,  J)  =  G_J  M  HJJ  ,  J  )-QH  )+0H 
’  <1,J)=HC1,  JJ/CFr. 
i .,=£?<♦!  bj  (nil  ,  J)  -OH) 

(Hll  ,  J)  ) 

LO  IS  I  =  2,M 


j)=iu*FM(i,j)*:i(i,j)»Tii^i,j)+c2i:,ji*T(i,j*i) 
1 ,  J) ♦ 34  (TVj  )  D/cS  n  ,77 

.0)  GC  TO  2 5 

sT5Tr»"JT7C  — 


■  C,  J)rO. 

"T?  TTTTjTTr 
GC  TO  15 


hii  *m  is  mr  wu.m  v 

w.v'U  -v,*,  «  ;,v  VO  IOC 


25  IF  (CH.lT.0.1  Hfl  ,j)=  (H  (I  ,jr-OHI*ORF*TiR - 

FF  =  0. 

IF  (J.LE.KF(D)  FF  =  i. 

IF  (NFIDH.EO.J)  FF-F  ( I  ) 

M*,  J)=HU  ,  J)/<1.-FF*<  l.-DF*)  > 

EN=EN*ftB3«H(I , J) -OH) 

ED=EO*A3S  (h(I,J) ) 

15  COM  I NUE 
OH  =  H (M  »  J) 

H(f,J)=(C*FH<M,J)+C2<MfJ)*-<M,Jtl)+C3<M,JJ*T(M-l,J) 
ltCilr.,  J) '7 (*S  J-l)  >/C5C  J) 

_ IF  (H(M, J) .LT.Q. )  GO  TC  ««5 _ _ _ 

rll”,  J)  =H(K,  J)*C5  (H,  J)/C 

T  <*,  J)  =  0. _ _ _ _ 


G <j  VO  35 

*•5  IF  COH.LT.O.)  H  Jn,  J)  =ORF*  J) -OH)  *OH 

T(t  , J) =H(H, J)/DFk 
EN=EN+  nBZ  j)-OH) 

ID=EO*/ c3(h(H,j)) 


35  COMIWUE 

5  COMIKUc 

OH=H(l,N) 

HC1,  N) = (O*  Fh  C1,K)*C1 (1 ,M  »T(2,N) +C4C 1, N) *T IlfN-l) 
l-Oi  *  SI  )  /C5  ( 1 , :«) 

H  ( 1  *  ■) )  =  0-  F  *  (H(l  ,n)-OH)  *0H 

1  (1  ,  .4)  =H  (1  ,N)/CFr. 

_  £(.-£'*♦>  fcS<H(  1,M-QH) _ 

El/=EC*-i.‘>£?<H(1iN)  ) 

:c  55  :=e,mi 

OH=H(I,U) 

_ Hi:,Ni--(:'tH(i,t,)^cni.n),T(iu,N)*C3c,Mii‘Ti!-i,n) _ 

< I »  N-l )  )  /C5  ( I  *  F! ) 

IF  (H(J,N)  .UT.0.  )  GO  TC  65 
HCI, 4)=H(I,N)*C5(I,N)/D 
T  C  ,  ;4)  =  0. 

IF  (  VF  (  N)  .  LT  .  1.  )  RF(N)  = l I-I r-h (I . N ) )  *DF«-1. 

GO  TC  55 

65  IF  (OH.LT.O.)  h  <  I ,  N)  =C  -.F  *  (  H  ( i  ,  M  -OH)  *OH 
5F  =  0. 

IF  <MF«I).GE.S)  FF=1 .  " 

T  (I  ,  '!)=H(I*N)/(l.-rr*  (  l.-OFR)  ) 

E*;=EN*  VET  (H  (I ,  ►.) -OH) 

”I-IC+  e:>(6(I|f\)) 

5  5  1?  •’  Tt.  'JF 

=  H  (  r. ,  H) 

H  C .  ,.>.)  =  (Li't  “«rl‘.T.)*:3  (M.kT*!  CH-l,;))+C4<h,N)  *7  (m,N-1)  )/C^Tm,n)- 
*  (.  ,;OsM(r.,N)/CrK 
75  !/0‘  T  IMJE 


iF  C  E»i.  I'  .EJ*E-  f.  )  GO  TO  3  5 


I  F  (  J.EQ.fF  ( i  )  ♦  1 )  HI  =F  (I  ) 
HJ  =  mMAX1  (0. ,H(I , J) ) 
S'i=Sl'i*--»uS«Dr<*<HI-HJ) 
Hl-HE  +  C  *  DS  *  Dr  *  (HI- H(I,  J)  ) 

PH(I,J)=T(I,J)/ST 

E'.-E‘;/EG 


•  r1  —  L>  •  *  j 

K  /3  .1415^2  7 
HE=FE*  2 . •NFN/3. 141532T 
f.E  =  SCF  r  (1.+S1) 


H-iTE  (2 ,  4C0)  7  F  ,r  E  ,  HE  ,  S  *  ,  i  T  ,  EF 
Ji  =  l  " 

J2=J1+ J° 

J2=«  MI>,0(  J?,N) 

wc  I T  E  (2,450)  (\F < J)  , J  =  J1 , J2) 
ro  5  J=J1,J2 
-S  ( J  )  =  (  J-.  5)  *CS*  EC 
W-I7E(?,500)  li-.F  < J), J= Tl  ,  J2 ) 
uC  11  J=J1,J? 


Nr  ( J )  = 

w-.n  e 

C3.*’  (1 
(  2  ,  c  5  0  ) 

lO  14 

J=J1 , J2 

'  F  ( J  )  = 

<CGf.*T  ( 

wt.  ITF 

(2,750) 

WE  *7  E 

<  2,5G3) 

DO  12 

I  -  1  ,  H 

r  =1.  ♦  <  i-iF-,5)  *Cr. 

WHI7F  (2,650)  i-»  (PH< I , J)  , J  =  J1 , J2 
J 1  =  J  2+  1 

IF  (Jl.LE.) >  GC  1(TT3 
IC  =  0 

IF  (  IS  .  LE\TT  IC^l 


rHII,  J) 

IF  (IS.L7.TT)  G C  TO  H)  - 

EG  K,I  (5E12.0) 

FC--M17  (TT5) 

EjrMT  <lhl,/,2(  1H0,/)  ,1H  ,34X,*HEA7  T  r  m  i\  S  F  E  r.  IN  PHASE  CH4VGE/T  HE) 


.ill  ENEKGV  STC  .E.GE  S  V  S7  E v  *  .  /  .  lh  0  x  .  *  Cl  S  F  N  J  H  d  F  7  *  .73.* [77:|."~ 


Wti  nil  I*  ■*** 


69 


•  r  i  "4 


1 X , * F I N  LEFFTH,rT=*,c'10.6) 

350  r  C  Mf  T  t///,lH  ,56X,*  w  <  T  F  i  <j  i.  F  i  OPEr.”I  E2*  » // ,  i  H  0 , 39X  ,  *P  PC  °E  *T  Y  • 

1  *  ?0X  ,  *  F  H  «*S  F  CH»  NGE  -l/TFr'It_*  ,  4  v,  «  CftiTFI  STT^  M* Tr  =  ItT *7/71^11". 2 7X  D 
l"uEN2I*  Y  ,LEH/n*«3",  30X,  ^  A.  3,  16X  ,F;>.  3,  /  ,  1«0, 27X,  *7  HE"  MAL  C  ONDUIT; ■ 

I  I .  v,  9TU/FT  /H-./F*  ,15X,F£. 5,15X,F8. 4,/,  lH0T2T5<T,lrPT?'  iiAl  "IFFUS7VITY 
lF-’--*2/H‘5",?5X,Ft.6,15XfFc.6,/,lH0,27X,*LATENT  H£^7  Or  FUSION, STU/i 
13N*,12X,Fk.3,///,lH  »52>.  ,mIJEr.*.7  I6K  COnTIOL  ?Al  i  f'.E  7  E  RS*  ,  777TH  ,5X 
l*NOCAl  ST).U3TUF:E«-*,I5,*  iOC  FS  ‘•fOIALt_Y*,I5,“  AKGULAk  MOOES*, 15,* 

II  ME  21  EPS*  ,  /  ,  1H0 ,5X,*nE?.ATICN  P"^rTF?TRr«  OVE^ -RELAX* T70 K  'VIZ' 
1Cc.=  «,F6.3,5X,»EF*0.-c  T0lEFAi2E=*,£12.4,*  MAXIMUM  ITER4T 10  NS =*  ,  1 5 


360  FO  .PAT  (///,1H  ,  5X,*  =  EFFc  c.,2£  TE  MPER*  U  ’c  ,  c  =  *  ,  F7T73 , 1 IX  ,  »EU51  CN/IN 
_ 1.lt:al  -EmFEkATUKE,F=‘,F?.3,12X,»S1£FaN  NUN3En.=*  ,  r12.  j) 


■QTHB 

FO-  Mil 

(  1  F  1 

,-T:Fi_=.,Fr..*.,3x,-cFF  P.f.iiur=-,r?.6,'3x,-HtAT  extracted: 

l*-,Fo.e 

,  3X,* 

SOLi  Z  VOuUFEF*  ,F  3.6,3  X  ,  *17  EFtA”!  ONS  =  *  ,  14, 3X, 

l*EKr.OR 

=-,F12.5) 

For m; t 

C  1H0 

,«  S/L  RADIUS*, 15F8. 4) 

500 

FC*m  hT 

(1H 

, * F I N  ►  NGLE* , 16F  3. 1 ) 

650 

FORMAT 

(IF 

,*Iu  WALL  TEMP* , 15F i, A) 

600 

FO:  MlT 

(IF 

1-TF) /T 

-*) 

650 

rO»  MAT 

(IF 

,4X,1  6r£.4) 

700 

^0.  h AT 

(IF 

,15Ffc.4) 

750 

rO- MAT 

(IF 

,*OC  WALL  TEMP*,15F8.4) 

5  00 

Fo  MAT 

(IF 

,10112.5) 

c  00 

FO-  T 

(  IF 

,1515) 

STOP 

ENr 

•  -.« *»•*»*  *^2,* 


M&y 
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UD121 , T200  , I030,CM60000,STCSA.P740579,BAND0U,KL121 ,229-2835 
COMMENT . ##********NO  DECK******* *** 

COMMENT, ***********93392*********** 

ATTACH, F, FIN, CY=1  . 

FTN, I=F  ,  L. »  R=3  »  B=FLGO  . 

RETURN, F. 

REWIND, FLGO. 

FLGO  ,  F'L  =  20000  . 


*EOR 

.035 

4.11 

181  . 

1310. 

350. 

. 2244971 

13.5 

501 .12 

8691 .133 

.7056356 

.  0625 

.  1195833 

.002604167 

.005416667 

.05447917 

1.93 

21  80 
200  8 
*EOR 
*EOF 

.00001 

30  6 

0  20 

Figure  D-2 

Control  Cards 
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